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ABSTRACT 

The existing disk instability model for dwarf nova eruptions is reviewed, 
in the light of recent progress in the understanding of angular momentum 
transport in accretion disks. It is proposed that the standard lower branch in 
the "S-curve" in the effective temperature-surface density plane may not exist. 
Rather, angular momentum transport may be suppressed in quiescence as a 
result of cooling. The model for superoutbursts is also examined, and it is 
pointed out that recent simulations strongly support the idea of a thermal-tidal 
instability. 
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1. Introduction 

The main goal of any discussion of limit cycles in accretion disks is to explain behaviors 
like dwarf nova eruptions, the outbursts of black-hole x-ray transients, and possibly the 
outbursts of FU Orionis stars. The first of these is best manifested in the century-long 
light curve of SS Cyg (e.g. Cannizzo & Mattei 1992; Fig. 1), while the second is extensively 
reviewed by Chen, Shrader & Livio (1997). At the next level, additional phenomena, like the 
"standstills" of Z Cam systems (e.g. AAVSO observations, Fig. 2), and the superoutbursts 
of SU UMa systems (Fig. 3) must be considered. 

It is generally believed that all of the above phenomena are related to disk instabilities, 
with the principal behavior (dwarf nova and x-ray transient eruptions) being explained in 
terms of a thermal-viscous disk instability model (see e.g. reviews by Cannizzo 1993; Osaki 
1996; Lasota & Hameury 1998). 

In the present work, I briefly examine some aspects of this and related models and I 
propose a few possible modifications. 

2. The Disk Instability Model 

2.1. Steady State 

In most of the studies computing limit cycles in standard, thin accretion disks, the 
conservation equations are written by averaging over the disk thickness. For dwarf nova 
systems these equations take the form (e.g. Pringle 1981): 
mass conservation 

2ttH£(-V r ) = M (1) 
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Fig. 1. — The light curve of SS Cyg on the basis of AAVSO observations (from Cannizzo & 
Mattei 1992). 
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Fig. 2. — The light curve of Z Cam on the basis of AAVSO observations (from Warner 1995). 
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Fig. 3. — The light curve of VW Hyi, exhibiting normal eruptions and superoutbursts (from 
Warner 1995). 
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conservation of angular momentum 
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conservation of energy 
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Here E is the disk surface density, V R is the radial velocity, M is the accretion rate, f2 
is the angular velocity, v is the kinematic viscosity and M*,R* are the mass and radius 
(respectively) of the central object. Combining equations (1) and (2) for a Keplerian flow 
gives 

m r /r.\5 

(4) 
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Thus, approximately we have M ~ z/£ ~ T e 4 jj, a relation that will become useful later. 
Equations (1), (3), and (4) are normally augmented by a description of viscosity, which in 
the dimensional-analysis "a prescription" (Shakura & Sunyaev 1973) is given by v = acsH, 
where cs is the speed of sound, H is the vertical disk half-thickness and a is a dimensionless 
parameter assumed to satisfy a £ 1 (see §2.5 and §3). 



2.2. The Basic Local Limit Cycle 

First, it is important to note that there are at least three fundamental timescales that 
are associated with standard accretion disks; dynamical, thermal and viscous. (A fourth 
timescale is associated with the propagation of transition fronts, see §2.4.) The dynamical 
time is the period of a Keplerian revolution. From hydrostatic equilibrium in the vertical 
direction, it is also the sound crossing time of the disk thickness (or the response time to 
a perturbation of vertical hydrostatics). The thermal timescale is the ratio of the thermal 
content to the local dissipation rate, and the viscous timescale is the time it takes material 
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to viscously drift inwards. These timescales are given by: 

, o-i R H 
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Here -D(-R) is the rate of viscous dissipation (per unit area) and all other symbols have 
their usual meaning. As we can see, since a £ 1 and in standard thin disks H/R <C 1, in 
such disks t dyn < t th < t vis . 

The steady state equations (1), (3), and (4) are normally complemented by equations of 
hydrostatic equilibrium, mass conservation, and radiative transfer in the vertical direction. 
The computations of the local vertical structure (e.g. Meyer & Meyer- Hofmeister 1981; 
Cannizzo, Ghosh & Wheeler 1982; Smak 1982; Faulkner, Lin & Papaloizou 1983) were 
found to produce a multi- valued thermal equilibrium function in the (Log T e /j; Log S) 
plane (see Fig. 4; from the discussion following eq. 4, T e ff can also be replaced by M or 
vY). Note that this represents a local solution, at a given radial distance in the disk. The 
upper branch of the "S curve" corresponds to a hot, ionized state of the gas, while the 
lower branch corresponds to a cool, neutral state. It is easy to see that the lower and upper 
branches are stable (e.g. an increase in the surface density, which results in an increase 
in the viscous energy production rate, leads to an increase in the effective temperature, 
and concomitantly in the energy loss rate), while the middle branch is not. Imagine now 
that the rate at which mass is being transferred from the companion star is such that it 
corresponds to a point on the middle, unstable branch for some annulus (see Fig. 4). In 
such a case, no local stable equilibrium is possible. Rather, because the rate at which 
mass is being supplied is higher than the rate at which it can be transported locally, the 
surface density will increase (along the lower curve) until the critical surface density S ma:E 
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Fig. 4. — The S-shaped thermal equilibrium curve in the surface density-effective 
temperature plane (adapted from Osaki 1996). 
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is reached, at which point the annulus will heat up on a thermal timescale, "jumping" to 
the upper branch. There, since mass is being transported faster than it is being supplied, 
the surface density will decrease (along the upper curve) until the critical density E miri is 
reached. At that point the annulus will jump back to the lower branch, thus completing the 
limit cycle. 

Before turning to time dependent evolution it is instructive to examine the physical 
reasons for the existence of such a multi- valued S-curve. 

2.3. The S-Curve 

The physical reason for the S-shaped equilibrium curve in the (z/E, E) (or (T e //,E)) 
plane is the partial ionization of hydrogen at T ~ 10 4 K (e.g. Mineshige & Osaki 1983; 
Pojmanski 1986; Cannizzo 1992; Liu & Meyer-Hofmeister 1997). 

This partial ionization produces two effects which give the S-curve its shape: (i) there 
is a steep dependence of the opacity on the temperature in the range T ~ 6000-10000 K 
(Fig. 5), which causes a change in the sign of dlog(i/E)/dlogE, and (ii) the adiabatic 
temperature gradient drops substantially (Fig. 6), thus driving convection. Both of these 
effects produce a local maximum in log E as a function of log T e ff (Pojmanski 1986; 
Cannizzo 1992; See Fig. 7). Here, however, arises a difficulty which in my opinion has 
not been fully appreciated by disk instability modellers. While the opacity peak is hardly 
significant (e.g. in Fig. 7 it appears only for rather high (a ~ 1), possibly unrealistic values 
of the viscosity parameter), it is the convection peak which determines the location, and 
to some extent the existence of the lower branch. This fact raises serious doubts about 
the reality of the lower branch. It should be realized that it appears almost certain now 
that angular momentum transport and energy dissipation in disks are governed by MHD 



- 11 - 




1.5- 
1.0 
0.5 




3.6 4.0 
log T (K) 



Fig. 6. — The mean molecular weight fi and the adiabatic index Va<2- Note the decrease in 
Vad due to partial ionization (from Cannizzo & Wheeler 1984). 
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Fig. 7. — The effective temperature vs. surface density for (from left to right) logo; = 0.5, 
0, -0.5, -1.0, -1.5 and -2.0. From Cannizzo (1992). 
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turbulence (e.g. Hawley, Gammie & Balbus 1996; Brandenburg, et al. 1995). Under these 
circumstances, and with the possibility that much of the energy dissipation occurs in a disk 
corona, the role of standard (mixing length) convection is questionable at best. I will return 
to this problem in §3. 



2.4. Time Dependent Evolution 

Many time- dependent simulations of the disk evolution have been carried out by a 
number of researchers; the most recent ones being by Cannizzo (1998), Hameury et al. 
(1998) and Menou, Hameury & Stehle (1998). 

For the time-dependent case, the conservation equations can be written as: 
mass conservation 

<9E Id 1 dM tr 

¥ = -^™ + 2^1T (6) 

angular momentum conservation 

d(Zr 2 n) i d _ i d ( 3 <m\ j s dM tr 1 

energy conservation 

c pS f = 2(Q+ - Q~) - 2»UrF r) - c„V r ^ - ^ d ±M (8) 
at r or or fx r or 

Here M tr is the rate at which mass is incorporated into the disk, j s is the specific 
angular momentum of the stream from the inner Lagrangian point, r tid is the tidal torque 
applied by the secondary star, Q + and Q~ are the rates of heating and cooling (per 
unit area) respectively, and all other symbols have their usual meaning. While some 
small differences in the exact expressions used for the various terms exist among different 
researchers, these do not appear to be crucial for the general behavior. 
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The broad-brush evolution can be described as follows. In quiescence matter piles up in 
the disk, until at some radius the surface density exceeds S maa; (Fig. 4). The corresponding 
annulus heats up on a thermal timescale. Since the viscosity increases with temperature 
(even for a constant a), this annulus starts to spread, thus initiating the propagation of 
heating fronts, which eventually bring the entire disk to a high state (the upper branch). A 
typical evolution of the surface density, the accretion rate onto the central object and the 
disk mass is shown in Fig. 8 (taken from Mineshinge & Osaki 1985 and Cannizzo, Wheeler 
& Polidan 1986). Two cases are shown; one in which the outburst begins at a large radius, 
and one in which it begins at a small radius. As can be seen in the figures, since most of 
the disk mass lies in the outer parts of the disk, outbursts starting from the outside produce 
a much steeper increase in the accretion rate. Since in the high state S(r) ~ r _1 , while 
^min ~ f i it is always the case that S first drops below S m j n at the outer disk radius. This 
initiates the propagation of a cooling front which transforms the entire disk to a low state 
configuration, thus terminating the outburst. In cases in which the cooling front cannot 
propagate (e.g. due to strong irradiation of the disk by the central source), the evolution is 
slower, because it has to proceed on the viscous timescale. This may produce flat-topped 
maxima in the light curve. 

The structure and physical properties of the transition fronts have been studied 
extensively both analytically and numerically (e.g. Meyer 1984; Papaloizou & Pringle 1985; 
Lin, Papaloizou & Faulkner 1985; Cannizzo, Chen & Livio 1995; Vishniac & Wheeler 1996; 
Vishniac 1997; Menou, Hameury & Stehle 1998). Broadly speaking, these works have 
shown that the speed of the fronts is of the order of Vp ~ cxhCs, where an is the viscosity 
parameter in the hot state (see below) and c$ is the speed of sound in the front (for a more 
detailed discussion see Vishniac 1997). The width of the front was found by Menou et al. 
(1998) to be proportional to the disk pressure scale height H. 
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Fig. 8. — The left two figures (8a) show how the disk responds to an outburst which begins 
at a large radius. The right two figures (8b) show an outburst starting at a small radius. 
The top figures show the evolution of the surface density, the bottom two figures show the 
accretion rate and the disk mass (adapted from Cannizzo 1993). 
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2.5. Viscosity 

Angular momentum transport in accretion disks is due to magnetically driven 
turbulence (e.g. Balbus & Hawley 1998; Godon & Livio 1998 and references therein). 
Consequently, expressing the anomalous viscosity by means of a fixed parameter does not 
represent adequately the physical situation. The following discussion should be viewed 
therefore as a time and space overaging of the process of angular momentum transport. 
Since the MHD simulations indicate that Maxwell stresses dominate over Reynolds 
stresses (e.g. Hawley et al. 1996), the viscosity parameter a is given by the appropriately 
overaged value of B R B^j 'Aiipc 2 S) where B R and are the radial and toroidal (respectively) 
components of the dynamo-generated magnetic field. 

From disk instability models of dwarf nova eruptions, it was generally found that in 
order to reproduce the observed amplitudes and timescales one needs to assume that the 
viscosity parameter in the hot state, an, is larger than the one in the cold state, ac, by a 
factor of 5-10 (e.g. Smak 1984). Otherwise, only relatively frequent, small oscillations in the 
luminosity are produced. The changes to a local S-curve, introduced by using an /etc ~ 5 
instead of a constant value, are shown in Fig. 9 (Hameury et al. 1998). 

One may ask, do we have any observational handle on the value of al 

From visual amplitudes and rates of decline (when compared to models) one obtains 
olh ~ 0.2 (Smak 1984). On the other hand, Livio & Spruit (1991) have shown that from 
the recurrence times of dwarf nova eruptions one obtains ac ~ 0.05 (see also Smak 1996). 

2.6. A Representative Model Calculation 

The most recent, extensive numerical models have been carried out by Hameury et al. 
(1998). Unlike in some of the previous calculations, these authors use an adaptive grid 
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technique and an implicit numerical scheme, which allows them to resolve the transition 
fronts. Hameury et al. also allow the outer disk radius to vary — a boundary condition 
which proves to be important for the obtained behavior (see below). The results of two 
typical calculations are shown in Figs. 10-11, for two values of the accretion rate. The 
main difference between the two calculations is the fact that for the higher mass transfer 
rate (M = 10 17 gs _1 ) the eruption is "outside-in" (namely, starts at a large radius with the 
transition front propagating inwards; these are also known as "type A" outbursts, Smak 
1987), while for the lower (M = 10 16 gs _1 ) it is "inside-out" (type B). This basic difference 
in the behavior is well understood. In the case of a high mass transfer rate, material 
accumulates in the outer disk faster than it drifts inwards, hence the surface density first 
exceeds the local T> max in the outer part. For relatively low accretion rates S ma2 . is first 
exceeded in the inner disk. The obtained behavior of the disk radius agrees well with 
observations (e.g. Smak 1984; O'Donoghue 1986). As the outburst is initiated and material 
starts to be accreted, it transports its angular momentum outwards, thus the outer radius 
expands (Livio & Verbunt 1988). 

As I noted above, Hameury et al. examined the effects of using as an outer boundary 
condition a disk radius, R ou t, that is allowed to vary (see also Ichikawa & Osaki 1992; Livio 
& Verbunt 1988), as opposed to a fixed R ou t- They found that when R out was kept fixed, 
the outbursts always tended to be inside-out, and they involved the accretion of a larger 
fraction of the disk mass. The reason for the latter behavior is easy to understand. When 
R out is fixed, the material that moves outward (with the angular momentum it absorbed) 
piles up at the outer edge increasing there the surface density. Consequently, all of that 
material needs to be accreted before the surface density drops below £ m i„ and the cooling 
wave is initiated. The light curves obtained under the two different boundary conditions 
are shown in Fig. 12. 
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Fig. 9. — The dashed S-curves correspond to constant a, a — 
The solid line represents a = an on the upper branch and a = 
Hameury et al. 1998). 



an (left) and a = ac (right). 
«c on the lower branch (from 
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Fig. 10.— The outbursts obtained for M WD = 0.6 M , a c = 0.04, a H = 0.2 and M = 10 17 
gs _1 (from Hameury et al. 1998). 
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Fig. 11 — The same as Fig. 10, for M = 10 16 gs" 1 . 
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Fig. 12. — The light curves that are obtained assuming a fixed outer disk radius (lower 
panel), and allowing the outer radius to vary (upper panel). From Hameury et al. (1998). 
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I should note that in addition to successfully reproducing light curves and the behavior 
of the disk radius, the disk instability model has provided a reasonable explanation for 
the observed distribution of systems in stability/instability regions in the (M, P or b) plane 
(where P or b is the binary orbital period; e.g. Smak 1989; Osaki 1996). In particular, the 
Z Cam objects which experience "standstills" (Fig. 2) are found to be right at the border 
between the unstable and stable regions. Thus, it is expected that small changes in M can 
carry them into and out of the instability region, which can account for the standstills. 

3. MHD Turbulence and Limit Cycles 

Since angular momentum transport in accretion disks is driven by MHD turbulence, it 
is appropriate to ask whether the thermal limit cycle model as described in the previous 
sections is consistent with such a dynamo-generated "viscosity." In particular, it is 
important to attempt to understand what could cause the existence of two values for the 
viscosity parameter, namely an and ac- In this respect, the following suggestion has been 
made (Armitage, Livio & Pringle 1996; Gammie & Menou 1998): the MHD turbulence is 
quenched in quiescence, leading to a low viscosity, while it is operative in outburst, leading 
to a high viscosity. 

The exact mechanism by which the MHD turbulence is suppressed is not clear at the 
moment, but there exist at least two possibilities: (i) when the disk cools, the magnetic 
field does not decay as fast (it decays on a superthermal timescale, e.g. Tout & Pringle 
1992). Consequently, the condition for magnetorotational instability (Balbus & Hawley 
1991), B\ < 24/irpcg, is violated, and hence the dynamo action is suppressed, (ii) When 
the disk cools, the magnetic Reynolds number Reu = csH/rj (where i] is the resistivity) 
can decrease below ~ 10 4 . Numerical simulations indicate that the MHD turbulence may 
be suppressed when Reu ~ 10 4 (Hawley, Gammie & Balbus 1996). 
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While I am slightly more inclined to favor the second possibility, the important thing 
for the present discussion is the fact that viscosity is decreased in quiescence as a result of 
cooling. 

If viscosity in quiescence is indeed low, this should have some obvious observational 
consequences: (i) little accretion should take place in quiescence, (ii) There should be little 
or no emission from the boundary layer and little UV emission from the inner disk, (iii) The 
temperature profile in quiescence should be very different from that of a standard disk (it 
should basically be flat), (iv) A UV delay (with respect to optical) should be observed in 
the rise to outburst. 

There are strong indications that all of these predictions are consistent with 
observations at least in some systems. For example, observations of all of the black hole 
transients show very low accretion in quiescence. Similarly, observations of the cataclysmic 
variables OY Car and Z Cha indicate that the accretion rate in the inner disk is lower by 
about a factor of 100 than that in the outer disk (e.g. Wood et al. 1986, 1989). Also, the 
temperature profile in OY Car is extremely flat in quiescence, while it is consistent with 
that of a standard disk in Z Cha in outburst. 

It therefore appears that there is some merit to the suggestion that the viscosity in 
quiescence is low. However, an important question arises, if material accumulates more or 
less in a ring (due to the low viscosity), how can inside-out (type B) outbursts be obtainedl 
Before I attempt to answer this question it is important to examine the observational 
evidence that inside-out outbursts do indeed exist. There are by now 15 dwarf novae with 
claimed positively identified outburst types (e.g. Smak 1996; Warner 1995, Table 3.6). 
Outside-in outbursts are characterized by a rapid rise, the existence of a UV delay; and they 
follow a broad loop in a color-color diagram. All of these properties follow directly from the 
fact that the outburst starts from the outer disk, where most of the mass is concentrated 
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(see §2.6). Examples for such outbursts can be seen in VW Hyi (Smak 1987; Fig. 13) and 
they are well reproduced in theoretical models. Inside-out outbursts are characterized by 
a slower rise, no UV delay, and a narrow loop in the color-color diagram (Fig. 13). All of 
these, in turn, are a consequence of the fact that the outburst starts in the inner part of 
the disk, where there is relatively little mass. Examples to this type of outburst are seen 
in SS Cyg and AH Her, and again they are adequately reproduced by models. One can 
therefore conclude with some certainty that outbursts of the two types do exist. 

I now return to the question of how can inside-out outbursts be reconciled with the 
suggestion of low viscosity at quiescence. In particular, I note the following, seemingly 
strange fact (e.g. Warner 1987): SU UMa systems show only outside-in outbursts while 
Z Cam systems show only inside-out outbursts, in spite of the fact that the accretion rate 
in Z Cam systems is more than ten times higher than in SU UMa systems! 

I would like to suggest the following possible solution. If we take as the critical (for 
suppression of MHD turbulence) magnetic Reynolds number Reffl ~ 10 4 , then it can be 
shown (e.g. Gammie & Menou 1998) that the critical temperature corresponding to this 
value (so that for T > T crit , Reu > R^m 11 ) 1S approximately given by 



It can easily be shown then that the condition T > T crit assumes the form £ £ const., with 



where P or b is the orbital period and f(q) is some function of the mass ratio of the two binary 
components. Now, it turns out that Cam ~ 2£,su UMa- Therefore, it is possible that the 
disk has a non-negligible viscosity in quiescence in Z Cam (and similar) systems and no 
viscosity in systems like SU UMa. This could resolve the puzzle of inside-out outbursts. 

I therefore propose as a possibility the following new scenario for dwarf nova eruptions 



T, 



74ook n°- 2 £' 



,0.2^0.038 



(9) 



crit 




(10) 



-24- 




i 



■ i t i ! > [ i I 
•0 2 0.2 04 6 

B-V 



35 

-< 

- 34 

xy 
k. 

33 

ft 

32 



MODEL 4 


I ■ ■■ 1 ■ 1 J ■ 

MODEL 9 


' % . 




















111' 


1 ■ i \ 



120 I2S 110 :JS 130 i55 160 165 



TIME MO* S) 



Fig. 13. — Examples of observations and model calculations showing a rapid rise and a broad 
loop in the color-color diagram for outside-in outbursts, and a slower rise and a narrow loop 
in the color-color diagram for inside-out outbursts (adopted from Smak 1987). 
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(see also Armitage, Livio, & Pringle 1996; Gammie & Menou 1998). In the disks of systems 
in which at quiescence Rtu <C Reffl (or some equivalent condition for the suppression of 
MHD turbulence), there is essentially no viscosity at quiescence. These systems experience 
outside-in outbursts. The initiation of these outbursts may be aided by the ring of 
accumulated material becoming unstable to the Papaloizou-Pringle (1994) instability (e.g. 
Rozyczka & Spruit 1993). 

In the disks of systems in which at quiescence Rem ~ R^m 1 ( or an equivalent condition 
allowing MHD turbulence), accretion takes place also in quiescence. In such systems the 
possibility for inside-out eruptions exists. The main implication of this new scenario is that 
a "standard" lower branch of the S- curve may not exist. This would be consistent with the 
doubts I raised about the existence of such a branch in §2.3. It is clear that the viscosity 
cannot be zero at quiescence in all systems, since some systems (e.g. HT Cas, WX Hyi) 
exhibit low states which are almost certainly associated with a reduction in M (these low 
states would not have been seen if there was no viscosity). 

4. Superhumps and Superoutbursts 

I would like now to discuss briefly another disk instability, the one producing 
superoutbursts and superhumps. Superhumps are periodic photometric humps that are 
observed in some dwarf nova systems, and they have periods that are longer by a few 
percent than the orbital period (Fig. 14; see e.g. Osaki 1996: Warner 1995 for reviews of 
their properties). Most of the systems that exhibit superhumps (the SU UMa systems) 
have orbital periods below the 2-3 hr period gap in the distribution of cataclysmic 
variables (TU Men is an exception). Some black hole soft x-ray transients also apparently 
exhibit superhumps (O'Donoghue & Charles 1996). In the SU UMa systems, superhumps 
are observed mostly during superoutbursts (see below), although in some systems (like 
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VI 159 Ori) they have been observed after superoutbursts (e.g. Patterson et al. 1995). 

Superoutbursts are eruptions that are brighter than normal by about 0.7 mag, 
and longer than normal by a factor of 5-10 (Fig. 3). The general impression is that 
superoutbursts start with normal outbursts (e.g. Vogt 1974; Warner 1995). 

The basic model for superhumps has been suggested by Whitehurst (1988), and further 
investigated by Hirose & Osaki (1990) and Lubow (1991). It involves a tidally driven 
eccentric instability of the accretion disk. The eccentric disk precesses at a period P pre c- 
The superhump period is identified as the synodic period between the precessing eccentric 
disk and the orbital period 

— = — —■ on 

PsH Porb Pprec 

Non-axisymmetric waves in the disk can be expressed in the form exp[i(k9 — £Qt)] 
(where 9 = <fi + fit, with being the azimuthal angle in the corotating frame). The modes 
are thus identified by the pair (k,£). Hirose & Osaki (1990) and in particular Lubow (1991) 
demonstrated the crucial role played by the 3:1 resonance in the disk in the excitation of 
the eccentric instability (see also Whitehurst & King 1991; Molnar & Kobulnicky 1992). 
Very briefly, a perturbation in the eccentricity (characterized by (1,0)) combines with the 
resonant effect of the 3:1 resonance in the tidal potential (corresponding to (3,3)) to produce 
a two-armed spiral density wave (characterized by (2,3)). This spiral wave ((2,3)), in turn, 
combines again with the tidal potential ((3,3)) to amplify the eccentricity ((1,0)). For the 
instability to grow, the 3:1 resonance radius must lie within the disk. This implies that 
this instability can occur only for mass ratios satisfying q = mijm\ £ 0.25. This explains 
the fact that superhumps are found primarily below the period gap or in black hole x-ray 
binaries. 

Surprisingly, it has proven much more difficult to find a convincing model for the 
superoutbursts. The history of the field is nicely summarized in Warner (1995) and Osaki 
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Fig. 14. — The development of superhumps in V436 Cen. From Semeniuk (1980). 
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(1996). The most promising model has been suggested by Osaki (1989). In this model, the 
following sequence of events is suggested to occur. As the accretion disk undergoes normal 
dwarf nova eruptions, in each outburst the mass that is actually accreted onto the white 
dwarf is smaller than the mass accumulated (via transfer) during the inter-outburst period. 
This agrees with the results of detailed simulations (e.g. Cannizzo 1993). Consequently, 
both the mass and the angular momentum of the disk are building up gradually, and the 
disk radius experiences a slow overall increase with each successive eruption (Fig. 15). 

At some point, the disk radius hits the 3:1 resonance and the tidally induced eccentric 
instability sets in, producing a precessing eccentric disk. Osaki assumed that the eruption 
that ensues at that point results in an enhanced rate of removal of angular momentum by 
the tidal interaction (compared to normal outbursts), so that a much larger fraction of the 
disk mass is accreted, thus producing a superoutburst. Until recently, all of the suggestions 
in Osaki's model have been verified by numerical simulations, except the assumption about 
an increased rate of removal of angular momentum. Most recently however, a simulation by 
Murray (1998) showed that this assumption is likely to be correct too. Murray simulated 
normal outbursts by an increased viscosity and showed that once the disk encounters the 3:1 
resonance it becomes tidally unstable and the tidal torques become much more efficient in 
removing angular momentum from the disk. He further showed that the resulting increased 
accretion rate is consistent with superoutbursts. 



5. Tentative Conclusions 

In recent years it appeared as if the question of dwarf nova eruptions has been fully 
solved, maybe with just a few details left to be sorted out. On the other hand, on the 
question of superoutbursts considerable skepticism existed in relation to Osaki's model. 
The enormous progress in the understanding of angular momentum transport in accretion 
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Fig. 15. — The time evolution of the accretion disk luminosity, radius, and mass (from top 
to bottom) in a simplified thermal-tidal instability model for superoutbursts (from Osaki 
1996). 
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disks on one hand, and the increased sophistication of numerical simulations on the other, 
have contributed to a renewed interest in both of these problems. 

On the basis of the material presented in the previous sections I am inclined to make 
the following tentative suggestions: 

(1) Dynamo generation of viscosity may be suppressed in some systems (and reduced 
in others) in quiescence, due to cooling. The difference in the angular momentum transport 
properties may be the main cause for dwarf nova eruptions and x-ray transient outbursts. 

A continuing serious effort will be required, in MHD simulations on one hand, and in 
observations of systems in quiescence on the other, to determine the viability of the above 
suggestion. 

(2) Superoutbursts are probably caused by a thermal-tidal instability, as has been 
proposed by Osaki (with the normal outbursts however occurring as in (1) above). 

The best observational test of this model can be provided by monitoring the behavior 
of the disk radius through a series of normal outburst and a superoutburst. 

This work has been supported by NASA Grant NAG5-6857. I thank the Isaac Newton 
Institute for Mathematical Sciences for its hospitality and John Cannizzo and Kristen 
Menou for helpful discussions. 
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